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We develop a theory of nonlinear cosmological perturbations on superhorizon scales for a single 
scalar field with a general kinetic term and a general form of the potential. We employ the ADM 
formalism and the spatial gradient expansion approach, characterised by ©(e™), where e = 1/{HL) 
is a small parameter representing the ratio of the Hubble radius to the characteristic length scale 
L of perturbations. We obtain the general solution for a full nonlinear version of the curvature 
perturbation valid up through second-order in e (m — 2). We find the solution satisfies a nonlinear 
second-order differential equation as an extension of the equation for the linear curvature perturba- 
' tion on the comoving hypersurface. Then we formulate a general method to match a perturbative 

'^'n ' solution accurate to n-th-order in perturbation inside the horizon to our nonlinear solution accurate 

, to second-order (m — 2) in the gradient expansion on scales slightly greater than the Hubble radius. 

The formalism developed in this paper allows us to calculate the superhorizon evolution of a pri- 
mordial non-Gaussianity beyond the so-called SN formalism or separate universe approach which is 
equivalent to leading order (m = 0) in the gradient expansion. In particular, it can deal with the 
■ case when there is a temporary violation of slow-roll conditions. As an application of our formalism, 

we consider Starobinsky's model, which is a single field model having a temporary non-slow-roll 
stage due to a sharp change in the potential slope. We find that a large non-Gaussianity can be 
. generated even on superhorizon scales due to this temporary suspension of slow-roll infiation. 
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PACS numbers: 98.80.-k, 98.90.Cq 



I. INTRODUCTION 



Recent observations of the cosmic niicrovifave background anisotropy show very good agreement of the observational 
data with the predictions of conventional, single-field slow-roll models of inflation, that is, adiabatic Gaussian random 
primordial fluctuations with an almost scale-invariant spectrum [U, Ql- Nevertheless, as the observational accuracy 
improves, it has become observationally feasible to detect a small non-Gaussianity in the data ^;-4]. In particular, the 
PLANCK satellite @ launched last year is expected to bring us much finer data and it is hoped that non-Gaussianity 
■ may actually be detected. As a consequence, non-Gaussianity from inflation has been a focus of much attention in 
I recent years 

To study possible origins of non-Gaussianity, one must go beyond the linear perturbation theory p^| - [l2j . The 
conventional models of inflation cannot explain an observationally detectable level of non-Gaussianity, since the 
magnitude of it is extremely small, suppressed by slow- roll parameters [31 . Then a variety of ways to generate a 
large non-Gaussianity have been proposed. They may be roughly classified into two; multi-field models that pro- 
^ , duce non-Gaussianity classically on superhorizon scales jl4l 4l9| . and non-canonical kinetic term models that produce 
■ - ■ non-Gaussianity quantum mechanically on subhorizon scales |20l422| . In particular, in the former case, the SN for- 
malism (23l - [25j turned out to be a powerful tool for the estimation of non- Gaussianity 0,[l3[lH. 

In order to parameterize the amount of non-Gaussianity of primordial perturbations, the nonlinear parameter /atl 
is commonly used 0] . This is related to the bispectrum of the curvature perturbation on ^ (23j , and is generally 
defined as 
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Here Vq and Bq arc the power spectrum and bispectrum of C, respectively, and they are defined in Fourier space by 

(Cfcia.) = (27r)3^3(fci+fc2)^7'c(fei), 

(CfciCfc^Cfc:,) = {2TTf5\k^+k2+k^)B^{k^MM). (1-2) 

Corresponding to the two different origins of non-Gaussianity mentioned above, the nonhnear parameter /jv^ can 
be mainly classified into two types; the local type, /jvL^'j which may arise from multi-scalar models on superhorizon 
scales, and equilateral type, /^^'', which arises from non-canonical kinetic term models on subhorizon scales.^ 
The local type is called so because it represents a local, point-wise non-Gaussianity given by 

a^)^CG{x) + \fsT\U^), (1.3) 

where C,g is the Gaussian random field. On the other hand, the equilateral form of the bispectrum is given by 

/^7'(fcl,fc2,fe3) = y^, (1.4) 

with the shape function Anl typically in the form [2l] |. 

■^NL - k^ik] + i^Y. ^^'^f ' ^i*^ if = fci + fc2 + fcs. (1.5) 

i i<j i^j 

Note that the sign convention of /nl here follows WMAP's sign convention and it is opposite to Maldacena [l3|. On 
the observational side, the current bounds on the parameter /at^ by WMAP seven years in are -10 < fj^l^^ < 74 
(95% C.L.) for the local form of the bispectrum and —214 < < ^66 (95% C.L.) for the equilateral form. By 

PLANCK 0, it is expected that non-Gaussianity of the level |/jvl|>5 can be detected 3. 

In this paper, we investigate another possible origin of non-Gaussianity, namely, non-Gaussianity due to a temporary 
non-slow roll stage on superhorizon scales. In order to investigate such a case, however, the SN formalism is not 
sufficient. The reason is as follows. On superhorizon scales, one can employ the spatial gradient expansion [13- 
[29| . which is characterised by the expansion parameter e = 1/{HL) representing the ratio of the Hubble horizon 
radius to the characteristic length scale L of the perturbation. The 5N formalism or the separate universe approach 
is equivalent to the leading order approximation, i.e., 0{e^) in the gradient expansion. It is valid at a slow-roll stage 
when local values of the inflaton field at each local point (averaged over each horizon-size region) determine the 
evolution of the universe at each point. In the context of perturbation theory, this implies one can ignore the decaying 
mode of the curvature perturbation. However, when the slow-roll condition is violated, the decaying mode cannot be 
neglected any more and the gradient expansion to O(e^) is known to play a crucial role already at the level of linear 
perturbation theory [sol - ls^ . Thus, to evaluate non-Gaussianity from a non-slow- roll stage of inflation, it is neces sary 
to develop a nonlinear theory of cosmological perturbations valid up through O(e^) in the gradient expansion [33H35|. 

This paper is org anized as follows. In Sec. [Hi we briefly review the spatial gradient expansion to O(e^) in the uniform 
Hubble slicing |33l435l |. In Sec. IIII[ we develop a theory of full nonlinear curvature perturbations on superhorizon 
scales valid up through O(e^) in the gradient expansion. In doing so, we introduce a variable that represents the 
nonlinear curvature perturbation as an extension of the linear comoving curvature perturbation. Then we derive 
an explicit expression for the nonlinear solution. Wc flnd a nonlinear second-order differential equation which the 
general solution satisfles. In Sec. IIVI we develop a general formalism for matching a n-th-order pcrturbative solution 
to the general nonlinear solution we found on superhorizon scales. Then in Sec.|V]we consider a special case when a 
linear perturbative solution is matched to the nonlinear solution. This applies to the case when the inflaton is still 
slow-rolling when the scale of interest crosses the Hubble horizon scale. As an application of our formalism, we study 
Starobinsky's model in Sec. IVII Sec. IVIIl is devoted to summary and discussion. 

II. GRADIENT EXPANSION IN UNIFORM HUBBLE SLICING 

In this section, we briefly review theory of nonlinear cosmological perturbations valid up through O(e^) in the spatial 
gradient expansion, following [35j . Throughout this paper we consider Einstein gravity and a minimally-coupled single 



^ A new type of f^L has been studied recently [26l |. called the orthogonal type. This may be generated from higher derivative terms in 
the action. 



3 



scalar field described by the action, 
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(2.1) 



where R is the Ricci scalar and X — ~g'^'^d^(j)d^(f>. The reason why the scalar field Lagrangian is denoted by P is 
that it plays the role of the pressure as shown immediately below. We assume is timelikc. Then the stress energy 
tensor of the scalar field may be put in the perfect fluid form, 

T^. = 2Pxd^(l)d,^ + Pg^, = ip + P)Uf,u, + Pg^, , (2.2) 

where the energy density p and the 4- velocity u^^ are given by 

p{X,4>)=2PxX-P, u^^-^. (2.3) 

y X 



Hereafter, the subscripts X and 4> represent derivative with respect to X and </>, respectively. The following relation 
among the first-order variations of P, p and (f) will be useful in the analysis below: 

5P = cl5p + pT5(j3, (2.4) 

where 

^^ = 2pJx^Px ^ r=.l(P,-c?,,). (2.5) 

Note that Cs is the speed of sound for a gauge- invariant scalar perturbation in linear theory [36j . 

We consider a theory of nonlinear cosmological perturbations on superhorizon scales developed in [13, H^- We 
employ the ADM formalism and the spatial gradient expansion in the uniform Hubble slicing. To make clear the 
relation between the standard perturbative expansion and the spatial gradient expansion, let us introduce the two 
numbers n and m associated with the two different expansions: The number n denotes the order in the standard 
perturbative expansion with respect to the amplitude of perturbation, and the number m denotes the order in the 
spatial gradient expansion, 0(e™). For example, linear cosmological perturbation theory corresponds to {n = l,m = 
cx)), and second-order perturbation theory to (n = 2,m = oo). The SN formalism corresponds to {n ~ oo,m. = 0). 
Our study corresponds to (n = oo, m = 2). 

In the ADM decomposition, the metric is expressed as 

ds^ = -a^dt^ + -fijidx' + P'dt){dx^ + P^dt), (2.6) 

where a is the lapse function, /3* is the shift vector and the Latin indices run over 1,2,3. We also need the extrinsic 
curvature Kij defined by 

= - ^ [dti,, ~ A/?j - D,(ii) , (2.7) 

where D is the covariant derivative with the spatial metric 7ij. Then the variations of the action with respect to 
a and /3* lead to constraint equations, namely, the Hamiltonian and momentum constraints, respectively, while the 
variation with respect to the spatial metric gives dynamical equations, which may be written as a set of first-order 
differential equations for and Kij. For convenience, we further decompose the spatial metric and the extrinsic 
curvature as 

7y = aV'Tii, 



K,^ = -K%j+A,^ , (2., 



1 

3' 

and deal with the first-order differential equations for and (K,Aij), together with the Hamiltonian and 

momentum constraint equations. 

In addition to the Einstein equations, we also have the field equations for the matter sector. We note that because 
we only have a single scalar field, the field equation is equivalent to the energy momentum conservation equation 
V^T^jy = in the present case. 
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In order to solve the Einstein equations, one has to fix the gauge condition. As for the choice of temporal gauge, 
we adopt the uniform Hubble slicing, 

K = -3Hit), H{t) = — . (2.9) 

a 

The spatial gauge condition will be specified later, by (|2.18p . 

Now we employ the spatial gradient expansion. In this approach we assume that the characteristic length scale over 
which the metric varies is much larger than the characteristic time scale over which the metric varies. To apply it to 
the cosmological situation, we introduce a flat FLRW universe {a{t), ^o(i)) as a background,^ and suppose that the 
characteristic length scale L of perturbations is longer the Hubble length scale 1/iJ of the background, i.e., HL ^ 1. 
Then we attach a small parameter e to each spatial derivative in the field equations and expand them in e. Physically, 
the parameter is equivalent to the ratio of the Hubble radius to the length scale of perturbations, e = 1/{HL). 

The background fiat FLRW universe (alt), poit)) satisfies the Friedmann equation and the equation of motion 

- ^Po, ^dt {a'Poxdt(l>o) - Po^ - 0, (2.10) 

where = SttGn, Po = piXo,^, Pox = Px(^o>o), Po4> = -P^(^o,0o), and Xq = {dt(j)of . Since the FLRW 
background must be recovered in the limit e — > 0, natural assumptions on the metric are 

v^^O{e), /3^ = 0(e), (2.11) 

and dt'jii = 0(e). Actually, for this last assumption, following the arguments in 0, [sj] , we assume a stronger 
condition, 

dt% - 0(e2). (2.12) 

This corresponds to assuming the absence of any decaying modes at leading order in the gradient expansion, namely, 
the absence of spatially homogeneous anisotropy. This is justified in most of the inflationary models in which the 
number of e-folds of inflation N is much larger than the number required to solve the horizon and flatness problem, 
N > 60.3 

Applying the above assumptions to the field equations, we can estimate the orders of magnitude of various quantities 
in the gradient expansion. In summary, including the assumptions, we obtain the following estimates: 

^ = 0(1), 7y=0(l), v^=0{e), /3^ = 0(e), 

X-O(e'), i„ -0(e'), 5-0(e'), - O(e'), P = 0{e'), 

9t% =0(e2), dti^^Oie'), + ^ 0{e^), (2.13) 
where x, S, p and ip represent the fluctuations in a, p, P and 4>, respectively, defined as 

X = a-1, 5=^^, p=P-Po, (^ = 0-00. (2.14) 

PO 

Note that, these fiuctuations may be non-vanishing at leading order in the gradient expansion in general. The 
advantage of the uniform Hubble slicing is that these fluctuations all become of O(e^) [sj. We also note that the 
form of p for the scalar held system is specified by the relation (|2.4p as 

p = Po{c%S + Taip) + 0{e^), (2.15) 

where c% = Pqx / {"^PoxxXq + Pox) and Fq = {P04, - c^soPo<p)/po- 

We now spell out the general solution valid up through O(e^) in the gradient expansion [s^. We attach the 
superscript (m) to a quantity of 0(e™). At leading order, the only non-trivial quantities are '0 and jtj, which are 
given by 

^ = L^'>\x'')+Oie^), (2.16) 



^ In this section, the subscript indicates the background quantities. 

^ Hamazaki [s^l solved the nonhnear equation on the leading order in the gradient expansion, under a more general condition, dt^ij = 
0((5c) = 0{e^), where 5c is a small parameter characterizing the amplitude of decaying modes. 
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and 



(2.17) 



where L^^\x'') is an arbitrary function of the spatial coordinates {x'^} [k — 1,2,3) and f\^\x'') is a (3 x 3)-niatrix 
function of the spatial coordinates with unit determinant, respectively. 

To obtain the solution at second-order, it is convenient to constrain the shift vector more strongly than indicated 
by ([21^ : We set 



/3' = 0(e3). 



(2.18) 



The choice of /3* = is called the time- orthogonal gauge. We note that setting /3* = 0{e^) (or even = exactly) 
does not fix the spatial coordinates completely [s^ . There remains 3 gauge degrees of freedom corresponding to the 
coordinate transformation, 



-> = P{x') ^ 0(1) , 



(2.19) 



If we fix this 0(1) part of the gauge, then the gauge is completely fixed at O(e^) accuracy. Note also that the spatial 
gauge condition (|2.18p with + /3* = 0{e^) in (|2.13p leads to the comoving threading condition: = at O(e^) 
accuracy. 

With the above choice of gauge, the general solution valid up to O(e^) was obtained in [35|. It is given by 



5 = 



6^2 (p( 



X 

lij 



-L_a(K.»./'„m + c<^>)+o,.»). 



6k2(po + Po)a2 



l + 3c 



ipo + Po)aJt, J ipo + Po)a 
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AO] , A2} 

Jij J ij 



2 K 



(2) 



^ Hit')xit')dt' ] + Oie^) 
* dt' /■*' 



a(t")dt" + 



(2) 



dt' 



+ 0{e% 



,(2) 



ait')dt' + Cl'^]+Oie*) 



where the dot (') denotes d/dt and 



1 



(^(0))4% 
1 



{my 



(2,("))4y(o) 



1 .(0)p 



(2.20) 



(2,(0))4j(0) 

2 ( 2d, In lnL(") - D^°^ d':"^ In^ 



(2a,lnLW9,lnLW -^f A(°)lnL(")) fjf> 



(2.21) 



are the Ricci tensor and the Ricci scalar of the Oth-order spatial metric ■ , 



(0) 



and = R [(i^"'')"'/''^)] is the Ricci scalar of the Oth-order spatial metric (L^°^)^/jj . Here note that the 'constants' 

were absorbed into and f-^, 
explicitly for later convenience. The choice of the initial time of integration, , will be discussed in the next section. 

/•(°\ 0(2) and Ci;: 



of integration, L'^^) and fl^ , were absorbed into L^^^ and flj , respectively, in [33, while in this paper we write them 

2hoii 

The 'constants' of integration m\ f^°\ O^^) and cff are not mutually independent due to the Hamiltonian and 
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momentum constraints. They must satisfy 



^(0) 

Jij 



AO) 
Si ' 



det(i^"^) 



"•'(O)^i 



1, 
0, 



•'(0)'-"ij 



(2.22) 



where f^^^^ is the inverse matrix of f^j and Z)'-'') is the covariant derivative with respect to /^^ 



(0) 



III. NONLINEAR CURVATURE PERTURBATION 



In this section, we define a new variable that is a nonlinear generalization of the comoving curvature perturbation 
up through O(e^) in the gradient expansion. We then construct it explicitly. We find that this newly defined variable 
satisfies a nonlinear second-order differential equation which is a generalization of the equation for the linear comoving 
curvature perturbation. In this and the following sections, we omit the subscript from the background quantities 
since there will be no danger of confusion. 



A. Assumptions and definitions 



As mentioned in the previous section, it is necessary to fix the spatial gauge to fix the metric completely. To 
do so, we assume that the contribution of gravitational waves to jij is negligible. In other words, we focus on the 
contribution arising from the scalar-type perturbations. Then at sufficiently late times of inflation, we may choose 
the gauge in which "fij approaches the flat metric, 

7ii (i -> oo) , (3.1) 

where in reality the limit t — > oo may be reasonably interpreted as an epoch close to the end of inflation. This 
condition completely kills the remain 3 gauge degrees of freedom up through O(e^) in the gradient expansion. It may 
be worth noting that one may include gravitational waves by relaxing the above condition to 

^^\n^)^J~>0 (t^oo). (3.2) 

Imposing this condition at all times is equivalent to the gauge chosen in [l3j . 

At leading order in the gradient expansion, one may define the nonlinear curvature perturbation by 

= e^? . (3.3) 

It is known that C on the uniform Hubble slices is equal to that on the comoving {— uniform scalar field) slices at 
leading order in the gradient expansion Q. However, at second-order in the gradient expansion. This equivalence 
between the uniform Hubble slicing and the comoving slicing breaks down. Furthermore, the very definition of the 
nonlinear curvature perturbation (|3.3p becomes inadequate as seen below. 

Let us derive the relation between (^h and to 0{e^), where and in what follows we use the subscripts H and c 
to denote the quantities evaluated on the uniform Hubble and comoving slices, respectively. We have obtained the 
solution (|2.20p in the uniform Hubble slicing (temporal), time- orthogonal (spatial) gauge, 

K^-3H{t), (3\t,x') = 0. (3.4) 

In general, we have to consider a nonlinear transformation between different time slices. However, since (pn ~ O(e^), 
the transformation to the comoving slicing, tpc = 0, happens to be just like a linear gauge transformation. The 
comoving curvature perturbation is obtained as 

Cc - Ch - + 0{e'). (3.5) 

<P 

Then the general solution for valid up to the O(e^) in the gradient expansion in the comoving slicing, time- orthogonal 
gauge, 

ipS,xl = l3cit,x') = 0, (3.6) 
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is written by the solutions of ip and (p in p.20p . whose explicit form will be given in the next subsection. Here note 
that remains the same at O(e^) accuracy under the change of time-slicing from the uniform Hubble slicing to the 
comoving slicing, 



(3.7) 



Now we turn to the problem of properly defining a nonlinear curvature perturbation to O(e^) accuracy. Let us 
denote the linear curvature perturbation on comoving slices by Tt^™. In the linear limit, the variable (c reduces to 
TZ^^'^ at leading order in the gradient expansion, but not at second-order. The comoving curvature perturbation in 
the linear limit is given by 



n 



Lin 



H, 



Lin 



where, following the notation in [38|, the spatial metric in the linear limit is expressed as 
with Y being scalar harmonics with eigenvalue fc^ satisfying 

(A + fc2)rfc = 0, 

and 



Y 



Y . 



These expressions in linear theory correspond to the metric components in our notation as 



2Hj^ Yij , 



(3.8) 

(3.9) 
(3.10) 

(3.11) 
(3.12) 



where we have used the definition of C given by (|3.3|) . In the above we have omitted the eigenvalue indices and the 
sum over the eigen modes for notational simplicity, but the existence of it is implicitly assumed. For example, Q = qY 
means 



Qit,x')=}_^qkit)Ykix'). 

k 

Thus to define a nonlinear generalization of the linear curvature perturbation 



of HlY and HtY . Our nonlinear C is an apparent natural generalization of Hj^^'^Y, 

HlY = C . 



(3.13) 

we need nonlinear generalizations 

(3.14) 



As for HtY, however, the generalization is non-trivial. Nevertheless, because it is the O(e^) part of 7^, the cor- 
respondence may be assumed to be linear. Hence, by introducing the inverse Laplacian operator A"-'^ on the flat 
background, we define the nonlinear generalization of HtY as 



HtY = E 



[aV~^5V^(ln7), 



(3.15) 



With these definitions of Hi,Y and HtY, we can define the nonlinear curvature perturbation valid up through O(e^) 
as 



(3.16) 



As clear from p.lSp . finding HtY generally requires a spatially non-local operation. However, as we shall see in 
the next subsection, in the comoving slicing, time-orthogonal gauge with the asymptotic condition on the spatial 
coordinates (j3.ip . we find it is possible to obtain the explicit expression for the nonlinear version of HtY from our 
solution (|2.20p without any non-local operation. 
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B. Explicit expression 

First we derive the explicit expression of C,c- Using p.Sp with p.3p . it is obtained from the general solution p.20p 



C, = + + fK{t) + fc{t) + 0{e% 
where = 21nL(o) and ^(2) ^ 2X^2) /^(o)^ The functions fxit) and /c(t) are defined by 



H 



pT4> 



iP + P> Jt, 



a{t")dt" dt' 



H 



2(p + P)a3 J,^ 



a{t')dt' 



HpTcj) 



rdt' 



H 



2(p+P)2a3 2{p + P)a^ 
The terms K^^'> and C(2) are spatial functions, and the former is the Ricci scalar of the Oth-order spatial metric 

/^(2)[£(0)J^ ^r(i{0))4/(0) 



(3.17) 

(3.18) 
(3.19) 



(L(o))5 

-2(2A£(")+<5*^9,£(o)a/W)e-2f<°) 



(3.20) 



while the latter C^^) is arbitrary for the moment. 

At leading order in the gradient expansion, £^^^ is the conserved comoving curvature perturbation, equivalent to 
the fluctuation in the number of e-folds 6N from some final uniform density (or comoving) hypersurface to the initial 
flat hypersurface (on which C = 0) at t = i^. 



(3.21) 



However, at second-order this is no longer the case, since Cc(i*) is not equal to Cc(oo) in general as clear from (|3.17l) . 
The 'constants' of integration ^^2) — 2L^^'' / L'-^'' and C^^) at 0{e^) characterize the total time variation of from 
t = to t = 00. 

Next we consider the explicit expression of in the comoving slicing, time-orthogonal gauge. As mentioned in 
the previous subsection, 7^ is the same for both this gauge and the original uniform Hubble, time-orthogonal gauge. 
Hence it is the one given in (j2.20p . 



where we have introduced the integrals. 



A{t) = 



dt' 



a{t")dt" , B{t) = 



* dt' 



(3.22) 



(3.23) 



The time-independent terms f^^^ and f^^'^^ are determined from the condition p.ip as 



Jij ~ "V ' 



=2F,^/M(oo) + 2C^fi?(oo 



n(2 



(2), 



F^f is given by ^!T2T^ with = 21nL(o), 



_2£(0) 



(2) 

and Cj satisfies the constraint equation. 



69^ ( 



(3.24) 



(3.25) 



(3.26) 



9 



Now, from the above expression of jij, we derive the exphcit expression of E' = HtY . The definition of E in p.lSp 
with ^(0) = 21nL(o) gives 



Hence we first need to evaluate the fohowing expressions: 



^3.(°)^) 



(2) 



The latter can be immediately evaluated from the constraint p.26p to be 



,3«W^(2) 



(3.27) 



(3.28) 



(3.29) 



As for the former, recalling that F^'^^ is the traceless part of the Ricci tensor of the metric e'^'^''"' f!^]'' , we can verify 



„2£<'" AO) 
lij 



Thus by substituting (|3^ into (|3^ with ((3?24l) . we obtain 



where we have defined 



= 3i/(2) + ^A(t) + ^i3(t) + 0(e4), 



i7(2) = __L (i^(2)A(oo) + (7(2)s(oo) 



Of course, the linear limit of Ec = HtY found above reduces consistently to H^™Y. Setting 



cif «iA-i(a.a,-i5,A)(7(^)^-i5(^)r, 



which is consistent with the constraint p.26p in the linear limit, and 



(2) 



we find 



with 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



Finally, we obtain the nonlinear curvature perturbation TZ^^ defined by p.l6p in the comoving slicing, time- 
orthogonal gauge as 



(3.37) 



where the functions fnit), fc{i)i A.{t) and B{t) have been defined in p.lSp . p.l9p and (|3.23p . This will be the basic 
variable to be matched to the solution in n-th-order perturbation theory (n > 1). The determination of f*-^-* and 
C''-^-' by this matching will be discussed in the next section. 
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f.-shift 



Here, we investigate the dependence of our general nonlinear solution TZ^^{t) on the initial time which appears 
in the integrals in the functions fK{t), fc{t): ^{t) and B{t). Apparently the solution (|3.37p depends on if v^'' and 
are t*-indcpendent. However, if wc match the pcrturbative solution whose initial condition has been fixed deep 
inside the horizon to our nonlinear solution on superhorizon scales TZ^^{t), it should not depend on the choice of i*. 
This implies that the 'constants' of integration (purely spatial functions) and C*'^' must depend on in such a 
way to cancel the t*-dcpcndence of the temporal functions fxit), fc{t), ^(t) and B{t). In particular, this implies 
that TZ^^ must be invariant under an infinitesimal shift of t^,, 



^* — > + (5t . 



(3.38) 



Thus the invariance of TZ^^ on the variation of t* gives consistency conditions that must be satisfied in the matching 
formulas derived in the next section. 

Let us consider the variation of TZ^^ with respect to this boundary time in the integrals, 



St^TZ^^ = 7^^L(^, + St) - 7^NL(^,), 

where 6t^ is the variation with respect to i,. From (|3.18p and p.l9p . wc find 

giit 



(3.39) 



StjKit) 



fc{t)a{'t* 



3^2 



5tJc{t)^-^-^5t + 0{5t% 



5t + 0{6t'), 



where we have defined 



9iit) = 



_ i/(l + 3c2) 



Similarly, we obtain from p.23p . 



2(p + P)a2 



(t), 92{t) 



HpT(j) 



2(p + P)2a3 



(3.40) 



(3.41) 



-a{t,)B{t)5t + 0{5f), 



(3.42) 



Hence the variation of TiJ^^ is given by 



fc{t)a{t^) - ^) /V 5t - ^(7(2) 5i 



3k2 



'a[t^)B{t)^^5t- 



1 (7(2 



3^2 



■5t + 0{5t^). 



(3.43) 



12 a^{U) 12 

Then requiring that TZ^^ be invariant fixes how £^2) ^ylA C''-2) (and hence _ff(2) j should transform under this variation: 

^{2) ^ C'(2)+if(2)a(i,)5t, 

L J 



H 



(2) 



(3.44) 



Before closing this subsection, we mention that not only TZ^^ but also each of Cc and Ec is i^-independent as well. 
This is a refiection of the fact that they are all gauge-invariant because the gauge has been completely fixed. 



D. Second-order differential equation 



In this subsection, we derive a nonlinear second-order differential equation that T^^'" satisfies at 0{e^) accuracy. 
For this purpose, we rewrite the integrals in the functions fK{t) and fc{t) defined in p.lSp and p.l9p in terms of the 
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quantity z commonly used in the literature |39j . 



We also introduce the conformal time 77, 



a 



p 



(3.45) 



d^=^, (3.46) 

and use t and 77 interchangeably. 

In Appendix [Xj using the background equations, we derive formulas that are used to change the forms of the 
functions fxit) and fc{t)- Using these formulas, we obtain 

where the subscript * indicates the quantity evaluated at t = (or i] = 77*). Further, in the analysis below, it is 
useful to adopt the same notation for the integrals in the above equations as the one used in linear theory by Leach 
et al. [l^l and re-express them as 



where 

Here = D{ri^), = F{'r]:t) and "H* denotes the conformal Hubble parameter H = dlna/drj at 77 = ry*. The 
functions ^(77) and -8(77) are the same as A{t) and B{t) defined in p.23p except that they are now implicitly assumed 
to be functions of the conformal time. Note that t ^ 00 corresponds to 77 in the conformal time. Thus the 
functions D and F vanish asymptotically at late times, D{0) = F{0) = 0. Here it is important to note that the 
function 1^(77) is the decaying mode in the long-wavelength limit (i.e., leading order in the gradient expansion) in 
linear theory, and ^(77) is the fc^ correction to the growing (i.e., constant) mode, 

D" + 2-D' ^0, F" + 2-F' + c2 = 0, (3.50) 

z z 

where the growing mode is assumed to have the form 1 + k'^F{r]) + O(fc^), and the prime denotes the conformal time 
derivative, ' = d/drj. 

Let us first recapitulate the solution C,c given by p.l7p 

Uri) = + + fK{il) K^'^ + fciv) + O(e') , (3.51) 

and Ec given by p.3ip . 

Ec = 3^(2) + ^A(77) + ^B{rj) + Oie") . (3.52) 

Adding these two and using the new expressions for fK{v) and /c(?7), we find that the functions A{r]) and -6(7/) cancel 
out to yield 

= iio}+ii2)^H{2)_^hF{rj)-F,]K^^^+[Dirj)-D,]C^'^+Oie''), (3.53) 
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where is given by (I3.32[) . and we have introduced the new spatial functions i^"^^ and C^^-* by 

(2) ^ H^C^^'I ~(2) 



(3.54) 



This is one of our main resuhs. The solution turns out to have a very simple form. In fact, as noted in the previous 
paragraph, the functions D{r]) and F{ri) have special meanings in linear theory: D{r]) is the decaying mode at leading 
order in the gradient expansion and F{r]) is the O(e^) correction to the growing mode, as shown in p.50p . This implies 
that, within the current accuracy of the gradient expansion, our solution TZ^^ satisfies the nonlinear second-order 
differential equation, 

where ^^'■^•'[X] is the Ricci scalar of the metric obtained by replacing with X in (|3.20p . In the linear limit, it 
reduces to the well-known equation for the curvature perturbation on comoving hypersurfaces, 

nY''" + 2^n^"'' - c2 A[nY"'] = o . (3.56) 

Equation p.55p may be regarded as the master equation for nonlinear superhorizon curvature perturbations in 
second-order in the gradient expansion. It should be, however, used with caution. For example, since it is derived 
under the assumption that the decaying mode is absent at leading order in the gradient expansion, a decaying mode 
solution obtained from the above equation with O(e^) corrections cannot be justified. Nevertheless, it would be 
interesting to investigate if a solution to p.55p with the right-hand side set to exactly zero can actually be a useful 
approximation to a full nonlinear solution on the Hubble horizon scales or even on scales somewhat smaller than the 
Hubble radius. 



IV. MATCHING CONDITION 



The general solution p.53|) for the nonlinear curvature perturbation on comoving slices TZ^^ has three arbitrary 
spatial functions, £^^\ and C'^^K Note, however, that the number of physical degrees of freedom are two since 
£^'^'> can be absorbed into i^'^K This is consistent with the fact that TZ^^ satisfies the second-order differential 
equation (j3.55p , or the fact that scalar- type perturbation has only a single field degree of freedom in the Lagrangian 
formalism. 

Physically these undetermined 'constants' of integration must be determined by the initial condition at a sufficiently 
early time when the scale of interest is well inside the Hubble horizon. There since the gradient expansion is not 
applicable at all, we have to resort to the standard perturbation theory. Let us assume that we have obtained the n-th- 
order perturbation solution under an appropriate initial condition. Let us denote this perturbative solution by 72.^'^''. 
Introducing a small expansion parameter S (not to be confused with the density perturbation) that characterizes the 
amplitude of perturbation, we may write 

^cxact(^^) = 72P"''('y) + 0((5"+i) , (4.1) 

where T^^'^^'^* is the exact solution. Our task is to match this perturbative solution to our nonlinear solution on 
superhorizon scales where the accuracy of the gradient expansion to second-order is sufficient. 

In this section, we perform this matching at t = or 77 = 77*. We denote the n-th-order perturbative solution by 
^pert_ choose the matching time such that the characteristic comoving scale of our interest k crossed the horizon 
about one expansion time before. That is, we are interested in wavenumbers that satisfy 

At and around the epoch 77 = 77* , we assume that both 72.^°''' and TZ^^ are reasonably accurate approximations to 
the exact solution. Namely, for some finite range of time interval around 77 = 77* , we have 

^ZP''\^^) = 7^:;"""'(77) + 0(5"+^) - 72^^(77) + 0(e^(5"+l) . (4.3) 

This implies that 72P"*(7/) satisfies the same second-order differential equation as (|3.55p up to the error of 0{e'^, 5"+^). 
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General formalism 



Since TZ^^ satisfies a second-order differential equation, the solution is completely determined once its value and 
the time-derivative arc given at a time. Thus, assuming we know the n-th-order perturbative solution, the matching 
condition at = ry, is given by 



NL\r 



^port 

in 



pert y I 



The first condition of (|4.4p leads to 
and the second condition of (|4.4p gives 

^(2) ^ jnr'YM 



(4.4) 



(4.5) 



(4.6) 



Note that the above equation means (TJ.p'^'^')' = O(e^). This is because we have assumed that there is no decaying 
mode at leading order, 0{e^), in the gradient expansion. 
Using the matching conditions (|4.5p and (j4.6p . we obtain 



3-H, 4 



(4.7) 



This is the general solution matched to TZ"^ at rj = rj^,. What we need to know is the final value of 7?,^ at sufficiently 
late times, r] Q {t ^ oo). ft is 



(4.8) 

where note the first equality which follows from the assumption p.ip : i?c —> at sufficiently late times. Parallel to the 
second-order differential equation p.55p for T^^'^, which is a natural extension of the well-known linear version (|3.56p . 
the above expression for the final value of Ti,^^ turns out to be a natural extension of the result obtained in linear 
theory in [30| . 



B. -independence 



In Sec. lIII Cl we considered the variation of HJ^^ under an infinitesimal shift of i», and derived consistency conditions 
on the variation of the undetermined spatial functions ^^^^ and C'^^-'. Here wc show that the final result (|4.8p obtained 
by the matching at t = is indeed i*-indcpcndcnt. 

If we take the variation of 7^^'^(0) given by (|4.8p with respect to the matching time t*, we obtain up to errors of 
0(e4,<5"+i), 



= 5r] 



(^perty 



.(-^pcrty/ 



pert \ / 



F' 



z' D, 



(^pcrty/ ^ 2-(7^^'*)' + -fK^^'inp"''] 



pert] 



iv* 



-2 if (2) [T^Pert] 



where 6r] — St/a^,. Thus for 7^^'^(0) to be t*-independcnt, 72,p°''' must satisfy 

.(^pcrty^f' 

z 



(7^P'='')" + 2^(7^^'')' + Ji^^''[7^P"*] = 0(e^^"+') at 77 = 77, 



(4.9) 



(4.10) 



This is exactly what we assumed for 72,^°''*: It should satisfy p.55p except for additional errors of 0(5"+^^). Hence we 
conclude that 7^^'"(0) is indeed -independent. 



14 



V. MATCHING LINEAR SOLUTION TO NONLINEAR SOLUTION 

In order to determine the nonlinear solution TZ^^ , we need to know the values of the n-tli-order perturbative solution 
T^-P""^* and its first time derivative at the matching time r] = rj^. In general this is a formidable task. However, if we 
consider the case in which the universe is in the conventional, slow-roll single-field inflation at the stage rj < rj^, the 
linear solution is a very good approximation on both subhorizon and superhorizon scales. Even if intrinsic nonlinearity 
is important for subhorizon scale quantum fluctuations, such as the case of DBI inflation, there may be cases in which 
the evolution near the horizon crossing time may be well approximated by linear theory. In this section, we focus on 
such a case, that is, the case in which the nonlinear solution TZ^^ and its first time derivative can be determined with 
sufficient accuracy by the linear solution T^.^;"'" at the horizon crossing time. 

Note that, while there is no problem in defining Fourier components 7?,^^ of the nonlinear curvature perturbation and 
thus the corresponding horizon crossing time r] = r]k, Fourier components with different k do not evolve independently. 
Hence, we should use the same matching time for all Fourier components. Otherwise, it would not be obvious whether 
matching conditions for different Fourier components are consistent with each other. This is the reason why we have 
introduced 77*. 

Thus the nonlinear solution should be obtained by the replacements, 

nr'M ^ n^c"\v*) + siM , nr^'M ^ n^^'M + .'*2(r;*) , (5.i) 

in the right-hand side of (|4.7I) . where si,2('7*) = 0{5^) are functions of rj^ and spatial coordinates. This boundary 
condition uniquely determines 7iJ^^ and thus its Fourier components ^^.^^(jy), provided that si_2(?7*) are specified. 
The nonlinear part Si^2('?*) of the matching condition is determined by requiring that the resulting 'R^\{'ii) and its 

time derivative 72,^^' (77) evolved backward in time do not include nonlinear terms at r/ ~ r/^. This requirement is 
nothing but a restatement of our assumption that the evolution near the horizon crossing time be well approximated 
by linear theory. 

Below we first briefly review the general linear solution on superhorizon scales obtained by Leach et al. [s^] ■ Then 
we spell out the nonlinear solution in terms of the linear solution. Finally, we derive the bispectrum from our solution 
assuming that the linear solution is a Gaussian random field. 



A. Linear solution valid up to 0{e^) 

In linear theory, the curvature perturbation on comoving hypersurfaces follows p.56p . As usual, we consider it in 
the Fourier space. 



7^^;^"" + 2-7e^;^' + cie = o . (5.2) 



. z 
z 

Real space expressions will be recovered by the replacement k"^ — > —A at the end of calculation. 

The above equation has two independent solutions; conventionally called a growing mode and a decaying mode. 
We assume that the growing mode is constant in time at leading order in the long-wavelength approximation or 
equivalently in the spatial gradient expansion. Then in terms of the growing mode solution, u, the decaying mode 
solution, u, can be given as (soj 

The general solution of a curvature perturbation is written in terms of their linear combinations as 

7^L-(77)=aI^-^.(77)+/3^'"«(?7), (5.4) 

where the coefficients a^™ and may be assumed to satisfy a^™ + (3^™ — 1 without loss of generality. Note that 
the assumption of the gradient expansion (|2.12p corresponds to the condition, 

= 1 - 4'" = 0{e) . (5.5) 



This means, as mentioned before, that the decaying mode at leading order in the gradient expansion has already 
decayed after horizon crossing. 
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First we solve for the growing mode solution. In accordance with the gradient expansion, we set 

CJO 



(5.6) 



n=0 



At leading order in the gradient expansion, the growing mode solution u^*^-* is just a constant. Then inserting the 



above expansion with ul^' =const. to the equation of motion (j5.2p gives iteratively 



-1 + 2-<+i 



As shown in [33|, Oik^) corrections to u"^^ can be written as 



,(2) 



,(0) 



.(2 



(5.7) 



(5.8) 



where the integrals Diif) and F(rf) have been given in p.49p . and c'^^ and are arbitrary constants. We fix the 
two arbitrary constants as Cp' = and = —k'^Ff,/D^ so that Wfc(77») = u^"-* holds at 0{k^) accuracy. Hence 



D 



■D{r^) + F(ri) 



(5.9) 



As for the decaying mode, since the coefficient /3j[;™ is already of 0{k^), we only need the leading order solution. 
Since we may replace D with D in (|5.3p . we immediately find 



(0)/ N l.^ 



Thus from (|5.9p and (|5.10p . the general linear solution valid up to O(e^) is obtained as 



Mri) 



D 



F, 



,(0) 



(5.10) 



(5.11) 



Note that TZ^^ir]*) = u^^'' while 7?.)^'"(0) = Q!fc'"u^°''. Thus if the factor \a^™\ is large, it represents an enhancement of 
the curvature perturbation on super horizon scales due the 0{e'^) effect, which happens when the slow- roll conditions 
are violated. This will be discussed in detail in Sec. VI. 

Here it is useful to consider an explicit expression for a^'" in terms of TZ^^j^ and its derivative at 7] = r/, . From the 
general solution given by (|5.4p . we have 



T^:k(:>l) = «^'"^fc(r?) + (l-4'")^;fc(ry), 



(5.12) 



where we replaced /3j[j™ by 1 — a^™. Using the definition ()5.3p of Vkii]), the above equations evaluated at 77 = 77* read 



We may solve these equations for ^. The result is 



ml 



(5.13) 



(5.14) 



J n=v* 



At 0{k^) accuracy, we have 



UkM uI°\t],) 



Oik') 



3n,F, 



Oik').. 



(5.15) 
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and hence 



C.fc 



+ k^F, + 0{k^) . (5.16) 

ri=ri. 



In order to relate our calculation with the standard formula for the curvature perturbation in linear theory, we 
introduce 77/^ (or tk) which denotes the time at which the comoving wavcnumber has crossed the Hubble horizon. 

Vk = -^; 0<r<l. (5.17) 

The power spectrum at the horizon crossing time is given by 

{nY:^{vk)KjM) = i2nfp^\k)s'ik + fe'), p^\k) = \nY:^M\' ■ (5.18) 

As shown in the Appendix |B1 we can show the final value of the linear curvature perturbation as 

7^c:fc(0) = = S^k^'T^cikiVk) + 0{k^), (5.19) 

where 

a^'" = 1 + a'^Dk - k^Fk (5.20) 



and 



1 TZ' 



,2/ 



J rik 



This explicitly shows that independent of up through 0{k'^). The formula ()5.19p will be used in the 

next subsection. 

The power spectrum at the final time is thus enhanced by the factor |a^™p as 

{n];:m'Jii%,m - {2nf\a]i^^fpj^\k)s\k + k'). (5.22) 

B. Matched nonlinear solution 

Using the linear solution of the curvature perturbation given by (|5.1ip . here we derive the nonlinear solution by 
matching the two at r/ = r/». The main purpose of the matching is to make it possible to analyze super-horizon 
nonlinear evolution valid up to the second-order in gradient expansion, starting from a solution in the linear theory. 
In particular, we would like to evaluate the bispectrum induced by the super-horizon nonlinear evolution. For this 
purpose, we need to have full control over terms up not only to O(e^) but also to 0((5^), where we suppose that the 
linear solution is of order 0{S). Therefore, the matching condition at r/ = 77* should be of the form 

n^'-'M = 7^L-'(^^)+,S2(77*)+0(e^5'), (5.23) 

where 

siM = 0{5^), S2M = 0{5^) (5.24) 

are functions of and spatial coordinates. While the linear solution TZ^™{ri) is considered as an input, i.e., initial 
condition, the additional terms, si(ry,) and S2 (??*), sltc to be determined by the following condition: 



17 



• The terms of order 0(5'^) in TZ^]^ and Tl^Y should vanish at r] = rj^, where 7?-^^ is the Fourier component of 
TZ^^ and r/k is the time shghtly after the wavenumber k = \k\ has crossed the horizon; —ki-jk = r <C 1. 

In other words. si(?7,) and S2(??*) represent the 0{S^) part of T^^'" and TZ^^' , respectively, generated during the period 
between the horizon crossing time and the matching time. 

The matching condition ()4.5p and (|4.6p is written explicitly as 

C^'' = -^(7^^^'"'(^^)+52(r^))+0(e^53). (5.25) 



The nonlinear solution (|4.7|) is now expressed as 



where 



The corresponding Fourier component is 



(5.26) 
(5.27) 



si.kiv*) 



- - [F. Fir,)] k'l^i^., 



S2,kM - [F. ~ F{tj)] 



+ Oie\6^). (5.28) 



Note that, as already stated, the Fourier components with different k do not evolve independently. 

By demanding that the terms of order 0(5^) in TZ^]^ and 7?.^^ should vanish at ?7 = i]k, si,fc(??*) and S2,k{v*) are 
determined as 



Vk 



z'{i')ci{r,")di^" + 0{e\5-') 



S2,k{V*) 



1 



z^(7/)cf(7,')A/ + 0(e^53) 



(5.29) 



Therefore, by substituting these to (|5.28p . we obtain 



where 



" dTj' 



Vk 



Using the linear solution of the curvature perturbation given by (j5.11|) . we have 



— (i-4i"-fc2^,)4°\ 



(5.30) 



(5.31) 



(5.32) 



Substituting these into ()5.30p . taking the limit 77 — > and using ()5.19p yield the nonlinear comoving curvature 
perturbation at the final time 77 = (or t = 00) given by 

1 

4^ 



7^c.k(0) = 5rX'fc"fe)-I-^A^(0)Af^[7^^''"fe)] + O(6^5^) 

J>fc(U)Afc 

This is the main result. The first term corresponds to the result of the 5N formalism, the second term is related to 
an enhancement on superhorizon scales in linear theory, and the last term is the nonlinear effect which may become 
important if J^fc(O) is large. 



= n]^M{Vk) - (1 - a]:"^)TZ]^Z{rik) 7-^.(0)^i'^[^c"(^.)] + 0{e\5^) 



(5.33) 
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Bispectrum 



111 this subsection, we calculate the bispectrum of our nonlinear curvature perturbation by assuming that Ti^^^irik) 
is a Gaussian random variable. We assume the leading order contribution to the bispectrum comes from the terms 
second order in The final result (j5.33p gives 



= G(fc)7^^;^"(77fe) + i^(fc) 



^^^{Ak'^ - 5.Xk''')T^Y%'MnY%.,{m'')S\-k + k' + k") 



where 



G{k) = = l + a'^Dk-k^Fk, 
H{k) = ij-fe(O) 



\Fk. 



(5.34) 



(5.35) 



This is independent of r|^, as should be. Note also that, in general G{k) may depend on the directions of k. However, 
in the present case we may assume the absence of such spatial anisotropy: G{k) = G{k). 

By assuming the Gaussian statistics for TZ^^]^{rik) with the power spectrum (|5.18p . it is easy to calculate the power 
spectrum and the bispectrum of ^. 

The bispectrum is expressed in terms of the Fourier transformation of the three point function as 



{aXk,Ck,)c = (2^)'Bc(ki, k2,k3) S^k, + k2 + ka) . 



(5.36) 



where {■ ■ ■)c means that it extracts out only connected graphs. With the help of (|5.34p . the three point correlation 
function of ( is at leading order calculated as 

{Ck^Ck,Ck,)c = (2^') X 

'G*{k,)G{k2)H{k;) {4(fc? + fcf) - 26,,k{kij 5(3)(ki + k2 + k3)|7^L|^(rfc)^7^L;^^(r,,.J|2 + perms] , (5.37) 

where a superscript star denotes a complex conjugate and 'perms' means terms with permutations among the three 
wavenumbers. The power spectrum of 72-J;''^(?7fc) is written as (j5.18p . Then we have 



Sc(fci,fc2,fc3) = G*{kr)G{k2)H{k:i) {b{kl + k^) - fc|} F.^"^(fci)P.^"^(fc2) + perms 
47r4 



Z-»3 7^3 Z-i3 
/i]^fi,2fi-3 



G*{ki)G{k2)H[ki) {^{kl + kl) - kl] fc|p4°^(fci)7'^"^(fc2) + perms 



where 



V'^\k) 



KkM 



We define the fc-depcndent /^v^ as 



Bc{ki,k2,k3) - lfNL{ki,k2,k3) [|ciLin5L'"|2p40)(fc,)p40)(fc2)+ perms 



247r4 



fNLiki,k2, ks) |ci^i"5fen'^4°^(^i)^r(^2)fc3' + perms 



(0), 



h.3 1.3 



If V^{k) does not depend on k then 



(5.38) 



(5.39) 



(5.40) 



fNL{ki,k2,k^) 



J2 \Gih)Gik,)\Hl 



G*{h) G(fc,)i/(fcfe){5(fc2 + fc2) _ kl]k 



(5.41) 
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VI. APPLICATION TO STAROBINSKY MODEL 

There are several known models in which the O(e^) effect is important. For example, a potential of the form, 



V = -M' 



l + B- 



(6.1) 



can lead to two separate stages of inflation with a temporary suspension of slow-roll inflation in between the two 
stages [13, EOl ■ The Coleman- Weinberg potential can also lead to the same feature as discussed in [4l| . A similar 
feature is found in a theory with a non-canonical Lagrangian [s^ , 



P = -Vq{1 + Vi(I)^)VT^ , with X ^ -g'"'d^,<l>d^(l). (6.2) 
In this section, we consider a simple model with a temporary non-slow-roll stage flrst discussed by Starobinsky , 



Vo + A+{(j) - M for0>0o, 
Vb -f- A_(0- 0o) forc/x^o. 



(6.3) 



where > A_ > is assumed. The advantage of this model is that it allows analytical treatment of linear 
perturbations as well as of the background evolution, provided that Vo dominates in the potential, \V — Vq\/Vo <C 1. 
If A+ ^ A_, and for (p initially large and positive, the slow- roll condition is violated right after falls below (j)^. To 
be specific, the field enters a transient stage at which </> w —3H(p until the slow- roll condition is recovered: 



■^fTA-J ""^+ for(/>>(/.o 
6Ho^-<^ -A_-(A+-A„)e-3ff»^* for(^<(^o 



(6.4) 



where At = t — to with to being the time at which (p = (j)Q and the Hubble parameter H is approximated by i?o = H{to) 
because |F - Vo| < Vq. Then for At > we have 



A+-A. 



Thus the slow-roll condition is violated during the stage 3HoAt < ln[(A+ — A^)/A^] if A+/A^ ^ 1. 



(6.5) 



A. Linear solution 

For Starobinsky's model, z defined by p.45p is given by 

A^e""^' 

ao 



for At <0. 



Substituting this into the integrals seen in (|5.20p . we obtain 

I i + r(^'' 



and 

























\ak) 



1 / 1 /ao 

6 V«fc 



ftp 



It^\t^\ for ?7fc < ryo (tfc < to) : 
5 6 ao I 



1 / 1 /ao 

6 V«fc 
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for i]k > 770 (tk >ta), 



(6.6) 



(6.7) 



(6.8) 
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where = l/(— iJo^fc) = k/{rH()), fco = flo^^o is the comoving wavenumber that crosses the horizon at i = to, and 
we have introduced 



T 



1 



(6.9) 



Thus if Aj^/A^ 3> 1, we have Dk ^ 1 and fcpFfc ^ 1 for wavenumbers k not too much different from kg. 

The perturbations with these wavenumbers must have left the horizon before the transition time rj — r/o- Then the 
vacuum mode function for the comoving curvature perturbation of Starobinsky's model before the transition is given 
by the standard formula for slow-roll inflation, 



where ^ = —A+/{3Ho). This gives 



and hence 
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for r]k < Vo- 

The amplification factor a^™ is given by (|5.20p with a^, Dk and Fk calculated above. Thus, for rjk < rjo, we have 



= 1 - 
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For T ^ 1 and /c/fco < r <^ min[l, (fc/fco)^^/'^] , the power spectrum of the asymptotic value of the curvature 
perturbation is 



(7^L;^(o)7^^|„(o)) = {2i,fPn5Hk + k'), 



Pn = |«r«(^fe)l 



2fc302 



2fc2 
5 fcp 



1 fc6 
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(6.14) 



Note that the inequality fc/fco < is required by rjk < rj^ < rjQ^ i.e., the condition that the linear solution be matched 
to the nonlinear solution no later than rip. This expression for the power spectrum is known to agree with the exact 
result very well even for Tfe^/fep 3> 1 [S^ as long as rjk < ?7o- 



B. Nonlinear solution 



In this subsection, we match the linear solution of Starobinsky's model to the nonlinear solution on superhorizon 
scales by applying the formulas obtained in Sec.|Vl We focus on the wavenumbers fc < fco- Recall (|5.34p which gives 
the final amplitude of the comoving curvature perturbation. 



a = G{k)nY:^{vk) + H{k) 



^|^(4fc'^ - 5.,k''k''')n];^,{i^k-)n];^.{m''W{-k + k' + k") 



(6.15) 



where G(fc) and H{k) are defined in ()5.35p . In the present case they are given by 



9 k'^ 



fc3 



5 kl 3(1 - ir) kl 



T 



r 3 — ir 
Q 1 — ir 



1 








5 2 










H 


5fc2 




V 6 


rkoj 


12 



(6.16) 
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As mentioned at the end of the previous subsection, the power spectrum given by ignoring the r-dependent terms in 
G(fc) agrees well with the exact result without using the long wavelength approximation. This strongly indicates that 
the r-dependence of G{k) is an artifact due to incomplete matching of the exact linear solution with an approximate 
longwavelength solution. This observation is supported by the fact that it disappears in the limit r — > 0. Similarly, 
the r-dependence of H{k) must be also an artifact due to incomplete matching of the linear and nonlinear solutions. 
In this case, however, the term fc/(rfco) does not disappear in the limit r — > 0. This is because there should be a 
sufficient lapse of time for the nonlinear solution to evolve before the transition time to erase the memory of small 
errors in the initial condition, implying that the accuracy of the nonlinear solution increases on larger scales k/ko <^ r. 
In short, we should take the limit k/ko ^ ?' <C min[l, (fc/fco)r^/'^] in (16.16^ to get rid of the ?--dcpendcncc to obtain 

Gik)^l~l^T^~T, H{k)^^T. (6.17) 



C. Bispectrum 

Finally, we estimate the bispectrum Bq and the corresponding non-Gaussianity parameter /ArL(A;i, A:2, fca) given by 
([08]) and (lOTj) . respectively. With the help of (|6T7l) . wc obtain 
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where, for T > 1 and k/ka < r <C min[l, (fc/fco)Ti/^] , we have 
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(6.18) 



(6.19) 



The fact that the power spectrum (|6.14p is in excellent agreement with the exact result [301 suggests that the above 
result (|6.19p based on the same approximation should also give a good estimate even for Tfc^/fcg ^ 1 as long as 
'7fc < Therefore Bq and Jnl may be evaluated as 
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(6.20) 



We expect the above expressions to be fairly good approximations for the entire wavenumbers in the range k/ko < 1 
and for any value of T as long as it is large, though it cannot be justified in the rigorous sense. 
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Setting fc = fci = fc2 = /cs, wc find 



1--—T 
5 kg 



9kl 



3.21) 



where f'^\^{k) denotes /^vl for an equilateral triangle. For a fixed T, it takes a maximum value (/^'^jmax = 54r/25 
at k'^/kQ = 5/(2r). Fig. 1 shows /^^^ as a function of y = \/Tk/ko for T = 10^. Note that f^l^ik) is always positive 
in this limit. 

For comparison, let us also take the limit of a squeezed triangle, k — ki — k2, fcs — 0. In this limit, /jvl becomes 



1--—T 

5 A:q 



9 fee 



1 - 7—?^ H r 



3.22) 



where f^l{k) denotes /nl for a squeezed triangle. Fig. 2 shows /^^^ as a function of y = VTk/kQ for T = 10^. Note 
that f^l{k) becomes negative for large T. 
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FIG. 1: /^i(fc) as a function of y = VTk/ko for T = 10^ 
To see the shape dependence of Bq, it is convenient to define a dimensionless function, 



1 

jj j4 fy>3 ^3 ^3 
y ■''l'''2'''3 . 



cyclic 



and identify the dimensionless variables as 

y = VTfci/fco, xi = l, X2 = k2/ki, X3 = /c3/fci, 

so that 

87r4(7'4"')'r3. 



Bf(fcl,fc2,fc3) 



5fcQ 



-;B;,(1,X2,X3). 



3.23) 



3.24) 



3.25) 



Without loss of generality, we can restrict our attention to the region 1 — 2:2 < X3 < 3:2 < 1. We plot By{l, X2tXz)x2x\ 
for y = 1 and y = 10 as a function of X2 and X3 in Fig. 3. For small y, we see that the bispectrum has a peak at the 
squeezed shape. On the other hand, for larger y, there is a positive peak at the equilateral shape as well as a negative 
peak at the squeezed shape. 
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FIG. 2: The same as Fig. 1 but for f^l{k). 




FIG. 3: By{l,X2,X3)x2x\ as a function of X2 and x-i for y = 1 (left) and y = 10 (right). The bispectrum has a peak at the 
squeezed shape {x2 = 1, X3 = 0) for y — 1, while it has a peak at the equilateral shape {x2 = 0:3 = 1) and a negative peak at 
the squeezed shape for y — 10. 



VII. SUMMARY AND DISCUSSION 



We have developed a theory of nonlinear cosmological perturbations on superhorizon scales for a single scalar field 
with a general kinetic term and a general form of the potential to second-order in the spatial gradient expansion. 
The solution to this order is necessary to evaluate correctly the final amplitude of the curvature perturbation for 
models of inflation with a temporary violation of the slow-roll condition. We have employed the ADM formalism 
and obtained the general solution for full nonlinear curvature perturbations valid up through second-order in the 
gradient expansion. We have introduced a reasonable variable that represents the nonlinear curvature perturbation 
on comoving slices TZ^^, which reduces to the comoving curvature perturbation TZ^"^ in the linear limit. Then we have 
found that TZ^^ satisfies a nonlinear second-order differential equation, p.55|) . as a natural extension of the linear 
second-order differential equation. 

Then we have formulated the matching of the nonlinear solution to a perturbative solution at n-th-ordcr, 7?,^°''', 
on superhorizon scales, and obtained a formula for the final value of the nonlinear curvature perturbation expressed 
in terms of 7?,^®''* and its time derivative at the time of matching. Since the evolution of superhorizon curvature 
perturbations is genuinely due to the O(e^) effect, our formulation can be used to calculate the primordial non- 
Gaussianity beyond the SN formalism which is equivalent to leading order in the gradient expansion. 
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Then we have considered the case when the hnear approximation is vahd up to the time of horizon crossing for 
wavenumbers of physical interest. In this case, we have developed a method to determine quantities corresponding to 
7?.P*^''* and its time derivative at the matching time in terms of the linear solution. 

As an example of such a case, we have investigated Starobinsky's model [42j in which there is a temporary non-slow- 
roll stage during inflation due to a sudden change of the potential slope. We have found that non-Gaussianity can 
become large if the parameter T « which characterises the ratio of the slope before and after the transition, 

is large. For T 3> 1, we have found that the non-Gaussianity parameter for the bispectrum fNL{ki, fc2, ^3) is peaked 
at the wavenumbers forming an equilateral triangle, k = ki = k2 = k^, denoted by f'^\^(k). It is found to be positive 
and takes the maximum value f'^\^{k) ~ 2T at y = \fTkjkQ ~ 1.5 where fco is the comoving wavenumber that crosses 
the horizon at the time when the potential slope changes. This implies that, even for a relatively small T, say for 
T = 10, it is possible to generate a fairly large non-Gaussianity /atl ^ 20 at wavenumber fc ~ O.Sfcg. 

Our formalism can be applied to many other interesting circumstances in which the slow-roll condition is temporarily 
violated. To mention a couple of examples, a non-slow-roll stage appears in a double inflation model ^l| or in a specific 
case of DBI inflation [3^ . It is of interest to investigate the non-Gaussianity in these models by applying our formalism. 

It is also of interest to investigate the case when there is a step in the inflaton potential instead of a change in the 
slope, which was proposed to explain the 'features' in the cosmic microwave background anisotropy [13, IS]. The case 
of time-varying sound speed for models with non-canonical kinetic terms (43 | may also deserve future study, since a 
rapid temporal variation of the sound velocity violates a certain type of the slow-roll condition. 

Finally, here we have focused on the case of a single scalar field. An immediate issue is to extend the present 
formalism to the case of a multi-component scalar field. We plan to work on this and hope to report the result in the 
near future. 
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In this Appendix, we derive some useful formulas that are used in Sec. IIIIDI to modify the apparent forms of the 
functions and defined in (PTS)) and ([XTOI . 

Using the background equations, 
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Appendix A: Useful formulas from background equations 




(Al) 



we obtain a useful formula by integrating the time derivative of (a/iJ), 




(A2) 



Then using the quantity z defined in p.45p 




(A3) 
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the above can be further transformed to 



dt \a{t) Jf.^ 
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where the subscript * indicates that the quantity is estimated at t = (or rj = rj^,). With the help of (|A1I) and (|A4p , 
we obtain 
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where the second term in the right hand side of the above equation can be rewritten as 
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and similarly, we obtain 
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Appendix B: t*-independence in Linear theory 

In the linear theory the curvature perturbation TZ^H^irj) in the Fourier space satisfies 



n]:f + 2-1^ + clk'n];% = o, 
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and is related to u^"'' up to 0{k^) as 
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In our paper we have assumed that 
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By inverting the relation (jB2[) and setting r] = rj^, u^*^'' is expressed as 
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where 
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- k^Fk + 0{k^) 



»)=»;. 



= <l- k^Fk 



^ Lin 



Vk 



,2(^^r)+o{k')}niX\vk), (B9) 



where 



F, - 



Vk 



^'iv') .nk 



(BIO) 



By using (jBip . it is easy to show that 
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This imphes that 
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Therefore, we obtain 



(B13) 



where 



5^- = 1 - k^Fk 
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and Fk is given by (jB10|) . This exphcitly shows that aj^^u'j^^ is independent of 77* up through 0{k^) 
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